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ABSTRACT 

We study the constraints that future multifrequency Cosmic Microwave Background (CMB) experiments will 
be able to set on the metal enrichment history of the Inter Galactic Medium at the epoch of reionisation. We 
forecast the signal to noise ratio for the detection of the signal introduced in the CMB by resonant scattering off 
metals at the end of the Dark Ages. We take into account systematics associated to inter-channel calibration, 
PSF reconstruction errors and innacurate foreground removal. We develop an algorithm to optimally extract the 
signal generated by metals during reionisation and to remove accurately the contamination due to the thermal 
Sunyaev-Zel'dovich effect. Although demanding levels of foreground characterisation and control of system- 
atics are required, they are very distinct from those encountered in HI-21cm studies and CMB polarization, and 
this fact encourages the study of resonant scattering off metals as an alternative way of conducting tomography 
of the reionisation epoch. An ACT-like experiment with optimistic assumtions on systematic effects, and look- 
ing at clean regions of the sky, can detect changes of 3-12% (95% c.l.) of the OIII abundance (with respect its 
solar value) in the redshift range z G [12,22], for reionization redshift Zre > 10. However, for Zre < 10, it can 
only set upper limits on Nil abundance increments of ^ 60% its solar value in the redshift range z € [5.5,9], 
(95% c.l.). These constraints assume that inter-channel calibration is accurate down to one part in 10"*, which 
constitutes the most critical technical requirement of this method, but still achievable with current technology. 

Subject headings: cosmic microwave background - intergalactic medium - abundances 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) radiation 
constitutes one of today's most useful tools in cosmology. 
The properties of its brightness temperature fluctuations carry 
information not only of the young Universe in which they 
arise, but also of the evolving Universe that CMB photons 
encounter in their way to the observer. Fir st and second 
gener ation CMB experiment s hke COBE jBennett et al.! 
11994ft. Tenerife ( Gutierre z et alJ l2000h. Boomerang 
(|de Bernai-dis et al.' 1999), Maxima ("Hananv et al. 2000) 
and WMAP (Scergelet al. 2003, 2006; Bennett etal. 2003; 
lHinshaw"et al. 2003, 2006), aim to characterize the proper- 
ties of the intrinsic temperature fluctuations, imprinted in 
the CMB d uring the epoch of Hydrog en recombination at 
z '^ 1,100, dSunvaev & ZeldovichllI970 ). Third generation 
CMB experiments, however, scan the small angular scales 
searching for secondary temperature fluctuations, introduced 
well after CMB photons decouple from matter during 
recombination. Among all secondary anisotropics, those 
induced by the thermal Sunyaev-Zel'dovich (tSZ) effect 
((Sunvaev & Zel'dovich ( 1980), caused by the distortion 
induced by hot electrons as they transfer energy to low energy 
CMB photons via Compton scattering) and the kinematic 
Sunyaev-Zel'dovich effect (( Sunvaev & Zeldovich 1972), 
caused by a Doppler kick induced by moving electrons via 
Thomson scattering) are of outmost importance. 

Indeed, the fact that the tSZ effect is a redshift-independent 
signal generated in hot dense electron clouds makes the tSZ 
a powerful tool to measure the abundance of galaxy clusters 
at all redshifts, which by itself can provide information about 
the presence and nature of Dark Energy in our Universe, 
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(e.g., iHaiman et alJ i2001l) '). T hird genera t ion ex periments 
like South Pole Telescope (SPT, "Ruhl et al." (2004")), and the 
Atacama Cosmology Telescope (ACT, Fowler etal. (2005)) 
will scan the sky at millimeter frequencies in an attempt 
to characterise the mass and redshift distribution of galaxy 
clusters in our Universe. 

Here, we exploit the exquisite sensitivity level and angular 
resolution of these forthcoming experiments to study the 
impact of metals on CMB temperature fluctuations. Metals 
are first produced by stars at the end of th e Dark Ages and 
during reionization (for a review see, e.g., 'Bark ana & Loebl 
(2001)), and, as shown in Basu et al. (2004), they provide 
some effectiv e optical de pth to CMB photons due to resonant 
transitions. iBasu et alJ ( i2004.) (hereafter BHMS) showed 
that neutral species such as C, O and Si, together with ions 
like CII, OIII or Nil should leave their imprint on the CMB 
temperature anisotropy by means of resonant scattering of 
CMB photons on some of their transitions. They found 
that the leading consequence of such scattering is a slight 
blurring of the original CMB anisotropics (by a factor T,y^^, 
with r^,^ the optical depth of a resonant transition) in the 
temperature maps and thus a contribution to the power spec- 
trum proportional to the primordial CMB power spectrum 
with constant of proportionality given by t^,^. Peculiar 
velocities of metals should introduce new anisotropics, which 
however should only be visible in the very large angular 
scales. BHMS demonstrated that, under ideal conditions, 
future multifrequency CMB experiments should be able to set 
strong constraints on the abundances of metals (lO"'*- 10"^ 
solar fraction) in the redshift range z G [1-50], opening a 
new exploration window on the end of the Dark Ages and the 
reionization history of our Universe. 

In this paper we build upon the work of BHMS paying 
particular attention to realistic aspects such as cross-channel 



calibration, point spread function (PSF) characterization 
and foreground contamination, aiming to assess their real 
impact on the sensitivity of future CMB experiments to metal 
enrichment history. We study what level of experimental 
systematic and knowledge of contaminants is required in 
order to achieve a given sensitivity threshold. Our approach 
is general, but we will focus on ACT's specifications (e.g., 
instrumental noise level, angular resolution and sky coverage, 
beam reconstruction etc.). Assuming a very accurate PSF 
characterization and inter-channel calibration for this exper- 
iment, we forecast ACT's sensitivity to the signal generated 
by metals and the constraints on reionization history that can 
be inferred from it. Current and future dust observations 
will help to correct for dust emission to significantly low 
levels in already relatively dust-clean sky areas. Also, 
because of the redshift-independent frequency dependence 
of the non-relativistic tSZ, an accurate subtraction this effect 
should be possible if more than two frequency channels 
are available. We find that the following specifications 
are required to extract the signal from a realistic survey: 
cross channel calibration error at or below the 10"'* level, 
PSF characterisation at 0.1-1% level, efficient foreground 
removal. Foreground and PSF characterisation requirements 
are realistically achievebable with an ACT-like experiment. 
The required accuracy in cross-channel calibration can 
certainly be achieved with Fourier Spectrometers, (Partridge 
[1995; Halpernetal. 1995). Such an experiment will be 
able to set strong constraints on the reionization history; 
in particular it should detect the presence of newly formed 
on at redshift z ^ 14 for metallicities of the inter-galactic 
medium above 3% the solar value, at the 95% confidence 
level. These constraints will provide a critical consistency 
check for estimates of the optical depth to reionization and 
constraints on reionization history obtained from CMB 
polarization measurements. In addition these finding will 
provide a direct test of our understanding of reionization and 
metal enrichment of the IGM. 

Since the measurement of the metallicity abundance is sen- 
sitive to a redshift window determined by the central fre- 
quency and the bandwidth of the detector, multifrequency ob- 
servations will provide a tomographic reconstruction of the 
epoch of re-ionization. This should motivate future experi- 
ments aiming to perform reionization tomography via CMB 
multifrequency observations. We offer guidance for these fu- 
ture experiments and quantify their technical requirements. In 
Section 2 we briefly review how the signal we are trying to 
measure arises and in Section 3 we build a realistic model of 
the total signal to be measured by experiments. In Sections 
4 and 5 we compare the signal generated at reionization with 
contaminants, and suggest different methods to extract it. In 
Section 6 we make predictions for an ACT-like experiment: 
we build a realistic foreground model and compute and dis- 
cuss the limits that can be imposed on the reionisation and 
enrichment history of the Universe. We conclude in Sec. 7 
discussing our results and prospects for future and dedicated 
experiments. 

2. CMB RESONANT SCATTERING ON METALS 

Following BHMS here we review how the signal of met- 
als arises. Let us assume that a CMB experiment is observ- 
ing at a given frequency iy„bs and that CMB photons inter- 
act with a species X via a resonant transition of resonant fre- 
quency f„. It is clear that the CMB experiment may be sensi- 



tive to such interaction only in the redshift range centered on 
Zrs = Vrs/vobs- 1, with its width depending on the frequency 
response of the detectors^. The resonant scattering of this 
species will generate an optical depth to CMB photons which 
can be computed to be (ISobolevll946h : 



Tx{z)=fr 



ire^ \rsnxiz) 
nieC H(z) 



(1) 



Here, nx is the number density of the X species, /„ is the 
absorption oscillator strength of the resonant transition, A„ 
its wavelength, H(z) is the Hubble function, e and m,, are the 
electron charge and mass, respectively, c is the speed of light. 
The effect of resonant scattering on CMB anisotropics can 
be more easily seen when considering a single Fourier A:-mode 
of the CMB temperature fluctuations. Such mode can be ex- 
pressed as an integral of various sources along the line of sight 
CHu & Sugivama 1994): 
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In this equation, written for the conformal Newtonian gauge, 
77 denotes conformal time, h.{rj) = f(7y)exp[-r(77)] is the visi- 
bility function and Tirj) = J* drj'f(rj') is the optical depth to 
CMB photons. For Thomson scattering, f = aajrie with a the 
scale factor, «« the electron number density and ctt- the Thom- 
son cross section. Arc accounts for the intrinsic temperature 
fluctuations, and vt, denotes the Fourier mode for baryon pe- 
culiar velocities; (p and if) are the Fourier coefficients of the 
scalar perturbations of the metric. 

As the CMB is being observed at a given frequency Vohs, the 
opacity due to the resonant transition can be approximated by 



TX=Tx Soirj-l-lrs), 



(3) 



where 77, 5 corresponds to z„, that is, the redshift at which ob- 
served CMB photons could interact with the species X via 
resonant scattering. When plugging this opacity into eq.Q, 
one obtains that original anisotropics have been blurred by a 
factor e~^^ , and that new anisotropics have been generated at 
the same redshift: 

In the limit of r^ -C 1 (low metal abundance), we can retain 
only linear terms in tx- Thus the change in the anisotropics 
(At-At„„^) reduces to a blurring factor {-txAt„„^) plus a 
new term proportional to v^, which is important only in the 
large scales (see Appendix A of BHMS for details). On the 
sphere, the blurring term is A{dT/Ti)) = -tx(ST /To)cmb and 
in multipole space, Aaim=-Txa™^ . Small scale experiments 
will be only sensitive to the blurring term, which will may 
be confused with other components present in the map. We 
address this issue in the next 3 sections. 

3. MODEL OF MEASURED TEMPERATURES 

The fluctuations measured at each pixel location by a given 
channel will be, in general, a sum of different components, 
namely intrinsic CMB fluctuations, tSZ-induced temperature 
fluctuations, metal-induced temperature fluctuations (which 

^ The frequency bandwidth of the detectors is assumed to be much larger 
than the Hne thermal broadening 



we aim to extract), galactic and extragalactic contaminants 
and instrumental noise. Hence, the signal received by channel 
j observing a frequency i/j is modeled as: 



T; = /> 



Wj* 



TcMB + T.sz.j + J2 Tm.> + E Tf .i + N,- 



(5) 

where M stands for sum over metals and F for sum over fore- 
grounds, N denotes instrumental noise, and all j indices refer 
to frequency vj; fj is an overall factor containing the calibra- 
tion factor of the given channel, and Wj represents its beam or 
point spread function (hereafter PSF). The symbol * denotes 
convolution. One can rewrite the Fourier counterpart of the 
previous equation as: 
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(6) 

This time we have split the peculiar frequency dependence of 
temperature fluctuations of tSZ ig,szj), foregrounds (grj) and 
metals (gMj) from the multipole {£) dependence (i.e. gxj = 1 
for a reference channel). However to keep track explicitly of 
the "real life effects", we will work with the ai's, defined as: 



4. POWER SPECTRUM OF DIFFERENCE MAP VERSUS 
DIFFERENCE OF POWER SPECTRA 

4. 1 . The difference map 

The difference map between channels j and / in multipole 
space is. 
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Since we assume accurate calibration and PSF character!- 
zation fjBl^ - 0,„^ ^ 2(ry,. - ^ + e>„ - e[j ^ 2(6 fj^ + Se^J 

and /,B^,„ gx.j-fiB^^ gx,i ^ gx.j-gxjjor any frequency- 
dependent component X. Hence, eq. MQ\ can be re-written 

as 
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In this equation, B^^ = wljiwlj"""^"^ ■■ 



(7) 
1 + e^^ expresses 



the departure of the real PSF from the model PSF, (i.e. if the 
PSF is perfectly determined, e^^^^ = and Bj^^ =1). The noise 

has been redefined as N/„, = i^ijiW/j"""'"'. Finally, if we set 
some frequency channel / to be the reference channel, we can 
write the calibration factor fj as 



fj = '^ + Vj, 



(8) 



where tjj denotes the error in the relative calibration with 
respect to channel k. As working hypotheses, we shall 
assume that both calibration and PSF characterization are 
very accurate, i.e., 77-' ^ 1 and je^^^J ^ 1 . 

The signal induced by metals can be rewritten as 



gMJ "f,„ - T^ff flfm , 
M 



(9) 



where the frequency dependence of the effective optical depth 
measures the amount of metals present in the redshift range 
probed by the observing frequency and the frequency of the 
resonant transition (rs), Zrs = t^rs/i^obs- 1- The lowest fre- 
quency channel will correspond to the highest redshift range, 
for which we will assume that the metal abundance is lowest, 
and therefore will be regarded as reference channel. We will 
attempt to measure increments of metal abundances with re- 
spect to such redshift by studying the difference map and the 
difference power spectra between different channels. 



with St^'cc ; 
eff 



eff eff 



This equation already shows the main limiting factor in our 
procedure: the accuracy of the cross-channel calibration. If 
Sfji is negligible, then it should be possible to filter the ^a^J„ 
to recover (5t''L and discriminate it from the remaining sig- 
nals, which should show, a priori, a different spatial pattern''. 
However, if Sfji is non zero, then it will give rise to a spurious 
signal of identical spatial properties to the signal we are after 
(-6T^'ff a™^), and hence indiscernible from it. Therefore, 
the calibration uncertainties set an unavoidable limit to 
the sensitivity of a CMB experiment that aims to measure 
the signal induced by scattering on metals during reionization. 



4.2. Difference of power spectra versus power spectrum of 
difference map. 

We next investigate the possibility of detecting the metal- 
induced signal in the power spectrum. There are two relevant 
quantities in this context: the difference of the power spectra 
computed at two different frequencies, SCj' = C;[Tj] -C/[T/] 
and the power spectrum of the difference map, C;[Ty-T,]. 

The latter consists on an average of \5a'^J\^, which, according 



to eq. JlU . reads 
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We are assuming that errors in the PSF characterization are imcorrelated 
with the intrinsic QVIB field, making the field ^^^lufil^,^ spatially indepen- 
dent of 4^" 
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(12) 
For third generation CMB experiments, instrumental noise 
will not be a limiting factor, and therefore here we have ne- 
glected the cross terms where they appear. Agi^ = gx.j- gx,i 
for any frequency dependent component X, and TZIY} denotes 
the real part of Y; F,, and F,, denote different foreground com- 
ponents. This equation simplifies considerably under the as- 
sumption that 5fji + 6ej'^j and St^'^c are much smaller than the 
typical amplitudes of tSZ and foregrounds. In this case, the 
average (J2m /(2^+ 1)) of eq. ([T2l becomes: 

(IKJ') « H',^\')(Agfszf + T.(\^L\')(Mf)' 
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(13) 



The last term in eq. ( I12> has been neglected because the 
correlation between tSZ and foregrounds becomes negligible 
when massive cluster are excised from the maps (see sec. |6}. 
That is, the power spectrum of the difference map provides a 
template for the dominant foregrounds in the case where the 
metal-induc ed si gnal is subdominant in Tj-T,-. We will show 
in Section (16. 2> that the spatial properties of foregrounds 
inferred from this statistic can be used to optimize the search 
for signatures of reionization. 

We next compute the difference of the power spectra com- 
puted at the channels considered, /' and /. Noting that 
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a straightforward calculation for |fl£„J^- Wi^V yields 
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The first term in the RHS of this equation reminds us 
of the importance of the calibration, since no signal will 
be measured if Sf^' ~ '^''"'ff The requirements on the am- 
plitude of the PSF error Se-'' are less strict than those of 
Sf-'', since its real part can flip sign and eventually will 
average to zero after combining different modes. Note 
that our signal is proportional to Sr^'^f, an d not to (St-''^)^, 
JHernandez-Monteagudo & Sunvaevl2005h . Since our analy- 
sis is restricted to small (I > 200) scales, the correlation be- 
tween the Integrated Sachs-Wolfe (ISW) component and the 
tSZ and extragalactic foregrounds can be neglected, because 
these signals are only relevant at large angles. Indeed, the 
tSZ and IS W correlation shows small (^ 3-1 uK) a mplitudes 
at I > 300, ("Hernand ez-Monteagudo & S unvaev'^005). The 
correlation between the tSZ and other foregrounds (e.g. ra- 
dio galaxies) is also negligible once bright galaxy clusters are 
excised out. If we further neglect noise (as done in eq. 12 
above), then we are left with 
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If tSZ is the only significant contaminant in the maps, then 
the £ shape of its power spectrum can be obtained from the 
difference map power spectrum, eq. JlU . and removed from 
l'^£mP~ I'^taP- Note that this is only possible because the 
spectral dependence of the tSZ is known to high accuracy and 
bright clusters have been excised. In more adverse scenarios, 
eg. din should always provide a rough estimate of the ampli- 
tude and shape of the foreground power spectrum present in 
\'^em\^~ I'^toiP- Note that the signal we are searching for has a 
definite well know ^-pattern, identical to that of C^*^*. 



5. APPLYING A TEMPLATE FIT 
Since we know what the spatial pattern of the signal we are 
trying to unveil is; Tj,[n] = -Sr^'f^T^'^^ln], with n denoting a 
position on the sky, and we are aiming to measure one single 
parameter (St^'^), it is justified to write the following model 
for the measured difference map: 



(T,. - T,)[n] = N[n] - <5T^^^T^^*[n] , 



(17) 



where N contains all the contaminating signals. 

If we assume that the signal N is gaussianly distributed"*, 
then an optimal estimator for St-''^^^ is given by dGorski et alJ 
Il99^ : 

/ry ry -.1 ^—IryCMB 

pre ji -1 _ '^^j ^'' »-' A n C\ 

iiLfTgff J - - .rj,CMByQ-ljCMB ' (^'^^ 

with C the covariance matrix of N (i.e. all the spurious signals 
regarded as noise), which is simply the Legendre transform of 

* If this is not tlie case, our matched filter of eq. llSt still provides an unbi- 
ased estimate of i5r^L 



eg. din . The formal error associated to E[(5t''L] is then given 
by 
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One must keep in mind, however, that the inversion 
of C might not be trivial, since its dimensions are 
Npx X Npx, with Npx the number of pixels of the map 
to be analyzed. If Npx is greater than a few thousands, 
the inversion of the covariance matrix poses numer- 
ical problems. We adopt the approach presented in 
iHernandez-Monteagudo. Genova-Santos & Atrio-Barandel^ 
( I2006.) of dividing the map into different equal-sized patches 
for which the covariance matrix inversion is doable, and then 
combining the estimates of 5t''^^ obtained from each patch. 
Therefore, for each patch /3, we obtain an estimate of the 
increment in optical depths St^'cc and an associated error 
o'.s^j''' given by 
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(21) 



We combine the Srj^ff estimates into a variance-weighted 
single estimate 
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whose uncertainty takes into account the fact that estimates of 
5t^'£ for different patches may not be independent: 
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The second term in brackets vanishes if different patches are 
uncorrected, i.e., C^''* = for /3i /^2- 

6. PROSPECTS FOR ACT 



The Atacama Cosmology Telescope (ACT, iFowler et alJ 
(yOOS)) will scan around 200 clean square degrees in the 
southern equatorial hemisphere in three frequency bands: 145 
GHz, 220 GHz and 265 GHz. Its angular resolution and sensi- 
tivity levels (FWHM ~ 1', with ^ 2/^K of instrumental noise 
per pixel) have been optimized to study the tSZ effect in clus- 
ters of galaxies. Using these technical specifications we make 
predictions on the sensitivity of this experiment to the pres- 
ence of metals at the end of the Dark Ages. We find that, be- 
yond the strict requirements of cross-channel calibration and 
PSF characterization, the main limiting factor will be the pres- 
ence of un-subtracted tSZ signal together with galactic dust 
emission. In what follows we will consider two cases, a pes- 
simistic and an optimistic case, chosen to span the reasonable 
range of possibilities. 



6.1. Foreground model 

At the frequencies of ACT, we consider three main fore- 
ground contaminants: galactic dust, tSZ and infrared (IR) 
galaxies. Radio galaxies should contribute with sub-dominant 
contamination at these relatively high frequencies, and their 
spatial properties should be close to those of IR galaxies. 

Dust: We build a model for dust emission based on the 
dust antenna te mperature map provided by the WMAP team 
jBennettet all 12003) and available at the LAMBDA site^; 
this model was produced using the data of ' Finkbeiner et all 
([1999). After converting it into thermodynamic temperature, 
we have computed its power spectrum outside foreground 
masks KpO and Kp2, (open triangles for Kp2, filled circles 
for KpO in Fig[^), and fit it with power laws. This should 
provide a fair estimate of dust emission at ACT's frequencies. 
In reality however, dust cleaning of the actual data would 
need to be done directly on the maps. The maps have been 
scaled to ACT's frequencies by using the intensity spectral 
index given in Bennett et al. (2003) (a = 2.2 ± 0.1). Note 
that this is a conservative choice: this spectral index overpre- 
dicts the dust emission i f compared to the sh allower index 
(a ~ 1.7) measured by Ponthieu et al.l (l2005h . We restrict 
our analyses to the cleanest third portion of the available 
sky in ACT's strip, for which the variance of dust-induced 
anisotropics is ^^ 4 % of the variance outside the KpO mask. 
Therefore, for our clean pixel subset we scale the dust power 
spectrum amplitude accordingly. We further assume that, 
after removing the existing dust model, the residuals vary 
from 10% (optimistic case) to 30% (pessimistic case) of the 
un-subtracted emission (in thermodynamic temperature). The 
dot-dashed line in Fig.([i];) shows the dust residuals in the 
power spectrum of the 265 GHz - 145 GHz difference map 
in the pessimistic case. 

tSZ signal: The first step consists in removing all clearly 
detected sources. In practice this means that all pixels show- 
ing tSZ signals above the ^ 4-cr level should be excised; this 
excludes practically all clusters above ^ IQ^'^Mq (see dotted 
line in Fig[0)). We are then left with the tSZ signal coming 
from smaller galaxy groups and the diffuse tSZ from a warm 
intergalactic gaseous component. When computing the power 
spectra of clusters and galaxy groups we consider both the 
Poisson term (C^, Atrio-Barandela & Miicket ( 1999)) and the 
correlation term (C^, iKomatsu & Kitavamal (Il999h ')^. These 
terms are given by 
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M„,. 



dM^[M,z]\yi{M,zt 



(24) 



and 



C^ 



(iz— - — P„{k = -— , z) X 
dz r(z) 

^dM^[M,z]b(M,z)y,(M,z) 
dm 



n2 



(25) 



For galaxy groups we chose M„„„ = 5 x IO'^Mq and M„,ax = 
5 X IO'-'M©; yi{M,z) is the Fourier transform of the tSZ an- 



^ LAMBDA URL site: \http://lambda. gsfc. nasa. gov I 

* The re-derivation of tiiese terms in a line of sight approach 
( Hernandez-Monteaeudo et al. 2006) yielded expressions of smaller ampli- 
tude. Since the tSZ is regarded as a foreground, we conservatively use the 
equations given above. 




Fig. 1 . — (a) Dust temperature power spectra as computed from the dust maps available at LAMBDA. Kp2 (empty triangles) and KpO (filled circles) masks 
have been applied, and the analytical fits used in our foreground forecast model, (b) tSZ signal present at 145 GHz. Once galaxy clusters above ACT's 4-cr 
threshold have been removed (dotted line), only galaxy groups and diffuse gas in superclusters and filaments remain. The total tSZ residual (filled circles) is 
dominated by the galaxy group contiibutions, both from the Poissonian (thin dashed line) and the coiTelation (thick dashed line). The contribution of diffuse gas 
is negligible (dot-dashed line), (c) Power spectrum of the 265 GHz - 145 GHz difference map assuming no metal signal (filled circles), and its three components: 
tSZ (dashed line), IR sources (dotted line) and galactic dust (dot-dashed line). 



gular profile of a cluster, 



y(n) = / liro-j-n, 



ksTe 



(26) 



where ar is the Thomson cross-section, kB the Boltzmann 
constant and Te, nie and Ue the electron temperature, mass and 
density, respectively; b(M,z) is the bias factor of lMo & White! 
( 11996;) relating the halo power spectrum to the underlying 
matter power spectra P„j(k,z)- Since we are interested in 
scales larger than the cluster, we took a gaussian profile for 
the gas distribution within these sources. Finally, we adopt 
a Press-Schechter mass function to characterise the mass 
and redshift distribution of collapsed objects. Since galaxy 
groups are smaller than galaxy clusters, the Poissonian 
term of the groups peaks at smaller angular scales than for 
clusters (see thin dashed line in Fig. ©))■ However, the 
distribution of galaxy groups is highly correlated, and this 
correlation dominates at large angles (thick dashed line in 
Fig. O?)). These two terms are responsible for virtually all 
the tSZ residual signal present in the maps after subtracting 
the brightest galaxy clusters (solid circles). We associated 
diffuse gas into massive (M > lO'^M©) structures that have 
turned around recently and hence still keep their linear size 
(superclusters). These structures are assumed to have an 
average overdensity of 6sc ^ 20 and a gas temperature of 0.2 
KeV. These parameters should over-predict the diffuse tSZ 
signal and hence provide an upper limit; still this is negligible 
compared to the emission generated in galaxy groups. In 
Fig©) the dot-dashed line gives the amplitude of the diffuse 
tSZ, which barely reaches the 1 /iK level at sub-degree scales. 

Infrared galaxies: Infrared galaxies present in ACT fields 
will be observed by BLAST^. Therefore, the major source 
of uncertainty should not be the amount of galaxies present 
within each ACT PSF, but how the measured flux at BLAST 
frequencies should be extrapolated at ACT frequencies. 
BLAST has three channels at 600 GHz, 860 GHz and 1200 
GHz, and depending on the galaxy's redshift these frequen- 
cies will probe the dust emission spectrum peak (high red- 

' BLAST URL shs:\http://chilel.physics.upenn.edu/blastpublic/index.shtml\ 



shift) or the long wavelength tail (low redshift). To estimate 
flux of these sources at ACT lower frequencies, we need an 
estimate specific intensity spectral index a in the long wave- 
length tail, (i.e. /^ oc ^" but a is expected to vary as a function 
of frequency). 

We thus build the following model for the infrared galaxies 
foreground signal and its residuals. We use the results of 
iKnox et alJ t.2004.) . which, based upon SCUBA** observa- 
tions, provides a spectral index of ascuBA = 2.6 ± 0.3 for the 
IR intensity over the range frequency 150-350 GHz. We use 
oiscuBA to extrapolate IR source amplitudes from SCUBA to 
ACT frequencies. We assign BLAST sources an effective 
spectral index of agiAST ~ 2 (shallower than that of SCUBA) 
to convert intensities between BLAST and ACT frequencies. 
This fixes our IR source model. We then introduce an error 
(SaBLAsr) in the spectral index of 5% (optimistic case) and 
20% (pessimistic case) of the model value of ublast, and 
regard the subsequent error in extrapolated intensity as and 
estimate of residuals IR source subtraction. We neglect any 
type of spatial correlation of IR sources and hence the IR 
source residuals show a Poisson spectrum (see dotted line in 
Fig.©, which corresponds to Sublast I ctBLAST = 0.2). 

We see that infrared galaxies contaminate ACT maps at 
small scales, particularly at 265 GHz. This level of contami- 
nation can be lowered by excising out the brightest galaxies, 
but here we prefer to compute the forecasts without mask- 
ing out more area. After all, at the angular scales at which 
the intrinsic CMB is dominant, the IR galaxy contribution is 
still below the tSZ residuals. Only at large angular scales the 
contribution from galactic dust is relevant. Fig.QJ shows our 
expectation for the power spectrum of the 265 GHz - 145 
GHz difference map (filled circles), assuming that there is no 
metal-induced signal, together with its three components for 
our pessimistic case. This is the most contaminated scenario, 
as in the 220 GHz - 145 GHz difference map the dust and IR 
emission is remarkably smaller, and the tSZ signal is reduced 
since the non-relativistic tSZ brightness fluctuations vanish at 

** SCUBA URL site: \http://www. if a. hawaii. edu/users/cowie/scuba/scuba. html\ 



218 GHz. 

6.2. One step further: removing the tSZ 

In Fig. QJ; we see that the tSZ signal due to galaxy 
groups and small clusters is the main contaminant. Since 
for most of these sources the electronic temperature is 
such that relativistic corrections are negligible, it should 
be possible to exploit the known frequency dependence of 
the tSZ temperature fluctuations in order to remove them 
from the difference map. This is feasible if the CMB exper- 
iment has more than two frequencies, as it is the case for ACT. 

For the sake of notation, let us refer to the 145 GHz, 220 
GHz and 265 GHz channels of ACT as the /, j and k bands, 
respectively. Let m and n denote any arbitrary pair of ACT 
bands. The / channel is probing higher redshifts, and therefore 
will be used as the reference channel, so metal-induced optical 
depths in the j and k channel can be written as r^-,^ = t' -,-,+ 

5t^'^^ and t*-,^ = r' r.^+(5r*r.^ respectively. Given any pair of 
ACT bands m and n, it is possible to build the following spatial 
template for tSZ signal: 



S.cz — 



^tsz 



A^; 



mn 
tSZ 



(27) 



rtSZ 



If the tSZ in band n is T„ = gtsz,n^tsz, then it is clear that 

^ mn 

S,52 will equal S,sz only if tSZ is the only component in the 
difference map, or dust and infrared emission, together with 
the metal-induced signature, are negligible. In such circum- 
stances, the spatial power spectrum of this template should 
be equal to that given in ea.( ll3> . and this step should provide 
a useful consistency check. Next we consider the difference 
map Ty-T,. For this difference map, it is possible to remove 

^ki " k j 

two different tSZ templates, namely S,^^ ^nd S,^^- That is: 



Note these estimates for /3/,' and /3,"" are computed over the 
same area on the sky, and the correlation between them arises 
due to the fact that foreground emission at different frequen- 
cies is correlated. This is also the cause for the correlation 
between (5R'"' and (5R^;, and for this reason we shall approxi- 
mate the correlation coefficient matrix of /3/,' and /3J"' by that 
of SW" and 5Rj ■. Finally, we must recall that we are not inter- 
ested on /3yJ or jSj!-, but on 5t''^ and St'^'^^. The latter are just a 
linear transformations of the former, and this makes straight- 
forward the transformation of the respective covariance ma- 
trices. It is clear that this procedure will be useful in case 
the tSZ component is far larger than the rest, since it removes 
most of the tSZ present in the difference maps at the expense 
of slightly increasing the level of the other contaminants. 



6.3. Sensitivity on Sr^'nr and constraints on IGM enrichment 
eff 

ACT will scan a strip of 2.5° width around the south 
celestial pole, at a constant declination of (5 ci -55.3°. The 
resulting ring covering 400 square degrees will intersect the 
galaxy, and hence the clean area will be smaller and located 
at high galacti c latitudes. In term s of dust emission, following 
the model of iFinkbeiner et aP tl999i) . we select the ^^ 140 
cleanest square degrees, and group them into 83 patches 
of 500 pixels each. This will make us sensitive to changes 
in the power spectrum in a multipole range I G [80,2000]. 
We use the HEALPix'' (Gorskietal. 2005) pixelization at 
resolution parameter Nside = 1024. The linear pixel size 
is close to 3.4 arcmins, and an effective convolution with 
a gaussian beam of FWHM= 4 arcmins is adopted for all 
components in all maps. We assume a PSF reconstruction 
uncertainty ranging from 0.2% (optimistic case) to 1% (pes- 
simistic case) (i.e., |(5e"™| ^ 0.002-0.01 VZ,m) and consider 
how the cross-channel calibration error affects the constraints. 



Rli=TJ-r-Agi-^i,. 



(28) 
(29) 



After introducing the dependencies on St-''^^ and ST'"ff, the last 
equations become 



K = 



K=- 



rCMB 



<5r^ff(l- 



^-eff(l- 



^^) + St'' 



^Sfsz 



^g'sz 



eff^ *./ 



tSZ 



) + Sr 



^slsz 



^^^8%, 



+ SRi'j,(30) 



+ <5R™.(31) 



The equations for R*' and Rj| constitute an equivalent system 
of equations (i.e. they no not add more information). The 
errors 5Ri'- and SR-'-: contain residuals of the tSZ subtraction, 

Kl Jl 

together with linear combinations of dust and IR galaxy sig- 
nals. From our foreground model, it is possible to estimate 
the covariance matrix of these errors. Further, if we apply a 
template fit introduced in the last part of Section (|5} on R™ 
and RJ"', we obtain estimates of the quantities 



/3/r 
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^-eff('- 



'5-eff(l- 
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)+K, 
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tSZ 
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(32) 



(33) 



When computing the CMB power spectrum we use the 
cosmological parameters that best fit WMAP's first year data 
JSpergel et alJl200l . and as a CMB template (X'^*'*) we use 
one sky realization from such power spectrum. Foreground 
contamination is introduced according to subsection l6.1l In 
what follows, in the pessimistic case we assume PSF recon- 
struction error of 1%, dust emission residuals of 30% in the 
maps (^ 10% in Q), 20% error in IR sources spectra index, 
no tSZ subtraction (after masking). In the optimistic case 
we assume \5('i",\ ~ 2 x 10~^, 10% error in dust amplitude 
modeling, SaBLAsr/ctBLAST, of 5%, tSZ subtraction down to a 
10% in the maps. 

In Fig.(|2^) we show the power spectra of the foreground 
residuals, for each of the two map differences, ji (or 265 
- 145 GHz, soHd lines) and ki (220 - 145 GHz, dashed 
lines). Thick lines correspond to the pessimistic case, thin 
lines to the optimistic case. We can see that the residuals 
power spectra become proportional to the CMB power 
spectra at / ^ 200, and that this fact is more visible when 
foreground subtraction is more accurate. This is due to the 
fact that PSF reconstruction errors Se^l are taken to have, on 



average, the same amplitude for all multipoles £'s. Hence 

I J J' CMl 



Jl „CMB\1 ^ (jCMB 



^ |http: //healpix. jpg.nasa.gov| 
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Fig. 2. — (a) Residual power spectra for the two map differences: 265 GHz - 145 GHz, (solid lines) and 220 GHz - 145 GHz (dashed lines). Thin lines 
con'espond to the optimistic case, thick lines to the pessimistic case (see text for details). The dotted hne shows the reference CMB power spectrum, (b) 
Minimum detectable St"'^ via a straight comparison of residual and CMB power spectra, i.e., ((5r"j?.)^ = |C^"/(2 C^"^^)], with C" the residual power as given 
in panel (a), c) Constraints that ACT will be able to set if inter-channel calibration issues are not important. In a scenario with early ionization powered by 
massive Pop III UV emission. Oxygen will be twice ionized and ACT will be able to measure increments of its abundance of ~ 3% its solar value between 
redshifts 22 and 14, and increments of ~ 12% its solar value between redshifts 22 and 12. (d) In a late reionization scenario in which Nil dominates the effective 
resonant scattering optical depth, ACT will detect increments in its abundance of ~ 60% its solar value between redshifts 9 and 5.5, and ~ 2.3 times its solar 
value between redshifts 9 and 4.5. 
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TABLE 1 

Sensitivity ON (5t^ 



Pessimistic case Optimistic case 
1.9 X Ur^ 2.2 X 10-* 



r. ki 2.7 X 10"^ 3.0 X 10-* 

e.ff 



We forecast at what precision the 5t'"S^ can be recovered 
via two different approaches. The first approach appHes 
the template fit outHned in Section |5] on the residual level 
given by our pessimistic of foregrounds, whereas the second 
provides the sensitivity on 5t"% when the template fit is 
applied on the optimistic case. The errors on 5t"% are 
shown in Table ^ for the two map differences considered: 
;■/ (220 - 145 GHz) and ki (265 - 145 GHz). Note that 
the ratio 5C'"" / (2C^'^^) can also provide a (less accurate) 
guess on 5t"'^^. Finally, let us remark that these three ways 
of estimating the increment of the metal induced optical 
depth, although not independent, can provide useful mutual 
consistency checks. For experiments scanning larger areas 
(i.e., SPT), the errors in the template fit method would 
roughly scale like the inverse square root of the area covered. 

The errors in Tableware given assuming that cross-channel 
calibration uncertainties are negligible. Thus if cross-channel 
calibration uncertainties {Sfji, Sfki) are worse than the values 
quoted in Table ^ then they will become the limiting factors 
in the measurements, (as already shown in eg. dlU ). Since 



abundances are proportional to optical depths (see eq.Q), 
it follows that errors in abundances are proportional to 
cross-channel calibration errors in case the latter dominate 
over all other sources of errors. These levels are reachable 
with current Fourier Spectrometer technology, JPartridge I 
fT995:.Halpern et al._1995-) . 

If indeed inter-channel calibration is not an issue, then by 
comparing Tabled with Table! in BHMS, we can infer that 
ACT will be able to set interesting constraints on the enrich- 
ment history. If Pop III stars ionized Oxygen twice, ACT will 
be able to detect increments of OIII abundances of ^ 3% its 
solar value between z ~ 22 and z ^ 14, and of ^ 12% its 
solar value between z ^ 22 and z '^ 12, at 95% confidence 
level, (see Fig.|2l;). Note that these redshift ranges are partic- 
ularly interesting, since they are close to the upper range for 
reionization redshift as indicated by WMAP (z,e ~ 1 1 ± 2.5 or 
5 < Zre < 15 at 2-0- level). Further, if reionization happens at 
the lower end of the WMAP values, and Nitrogen dominates 
the effective optical depth to CMB photons, then ACT will set 
constraints on the abundance of Nil, since it can measure in- 
crements of ~ 60% its solar value between z ~ 9 and z ~ 5.5, 
and of ^ 2.3 times its solar value between z '^ 9 and z ^ 4.5, 
(95%c.l., seeFig.Ell). 

7. DISCUSSION AND CONCLUSIONS 

The probability that a CMB photons is scattered by a given 
metal species depends primarily on two things: how abundant 
such species is and how likely the transition is. Indeed, as 
already shown in eq.Q, the Sobolev optical depth for the 
resonant transitions is proportional to the number density of 
(niiz)) the intervening species and the so-called oscillation 
strength of the transition (/,). This means that for fixed metal 



abundance, the oscillator strength, the transition wavelenght, 
and the relative population of the lower level involved in the 
transition will rule which lines are more visible. From Table 
1 in BHMS, the transitions giving rise to higher optical depth 
con-espond to OIII (88.36 ^m, 51.81 ^im), OI (63.18 ^im), 
CII (157.74 /im) and Nil (205.30 /im). For the frequency 
range of ground based experiments, O and N are the elements 
more likely to generate most of the optical depth to CMB 
photons. (Note that, although CII shows a strong transition, 
current Pop-Ill yield models predict low C/N abundance 
ratios, ( Mevnet. Eksti'om. & Maeder 2006) ). 

In the presence of only one species, say O, the technique 
presented here measures the increase of O abundance 
between some higher redshift bin (corresponding to the 
lowest frequency of the experiment) and lower redshift 
bins (corresponding to lower frequency bands). Under the 
additional assumption that the metal abundance at the highest 
redshift bin is negligible, the method yields constraint on 
the species abundance as a function redshift. In presence 
of different species (say O and N), this technique will be 
sensitive to the integrated effect of the different species. 
Note however that the same frequency band probes different 
redshift for resonant transitions of different elements. For 
example for OIII 145 GHZ corresponds to z ~ 22, 220 GHz 
to z 2± 14 and 265 GHz to z ~ 12, but for Nil, 145 GHZ cor- 
responds to z ~ 9, 220 GHz to z ~ 5.5 and 265 GHz to z ~ 4.5. 

To interpret these measurements in light of the reionization 
process, we need a model of metal production for pop-Ill 
stars and an estimate of star formation rates at these redshifts. 
There have been rece nt advances in the computation of the 
yields for pop III stars iMevnet. Ekstrom. & Maedeil2006ft : a 
more careful modeling including the effect of rotation in the 
evolution of massive, low metallicity stars makes predictions 
for the properties of low metallicity halo stars in better agree- 
ment with observations. Star formation rate estimates are 
somehow dependent on the initial mass function, however the 
yields are relatively insensitive ( Mevnet, Ekstrom, & Maeder 
120061) to the mass of the star for massive (M > IOMq) stars 
and pop-Ill are expected to be massive (Abeletal. 2003)- 
Therefore a constraint on the metal abundance obtained from 
measurements of the optical depth for resonant scattering, 
will translate into a constraint on the number of pop-Ill stars. 
Different reionization scenarios produce a different amount 
of pop-Ill stars as a function of redshift. While a detailed 
calculation along these lines will be presented elsewhere, 
we can show the potential impact of this technique and its 
complementarity with other probes of reionization with the 
following example. The nature of ionizing sources is not well 
understood (e.g. Barkana & Loeb (2001)). A detection of ^ 
12% solar metallicity of oxygen at z ^ 12 would imply that 
pop-Ill stars played a major role in reionization and thus rule 
out the scenario of reionization by mini quasars. 

In addition, the method presented here of tracing the reion- 
ization history with a measurement if the IGM enrichment, is 
complementary to other probes of reionization. For example 
the CMB is sensitive to ionizing photons produced during 
reionization and their escape fraction, and the horizon scale 
at reionization redshift. H 21 cm radiation probes neutral 
Hydrogen and the Ly-a forest has provided measurements of 
metallicities at z ^ 4-5. All these methods rely on different 
physics and thus provide a invaluable tool for understanding 



the physics of reionization. In addition they are affected 
by foregrounds and systematics of very different nature, so 
their comparison provides a useful consistency check for the 
results. 

The forecasts presented here are for an experiment de- 
signed with another goal in mind. It is interesting to ask how 
well an ideal experiment or a dedicated experiment could 
perform. A dedicated experiment would observe through 
the same PSF at different frequencies, in order to minimize 
the effects of PSF reconstruction errors. An self-calibrator 
device would assure accurate inter-channel calibration. 
The set of observing frequencies should be a compromise 
between accurate tSZ subtraction (this forces the ratios 
giszj/gtsz.j, gtsz,i/gtsz.k to be /or from unity) and the choice 
of an interesting redshift range in terms of cosmic enrichment 
history (which is determined by transition frequencies of the 
metal species dominating the resonant scattering). Detector 
spectral widths verifying /S.v/v J^ 0.1 should assure that 
resonant scattering is dominant over collisional emission (see 
discussion in BHMS). 

In our optimistic model for foregrounds, if the PSF 
reconstruction errors are smaller or equal than 5 x 10"^, 
then foreground residuals are the limiting factors and the 
constraints on the abundances improve by a factor of three 
compared to the case considered here, for the same sky 
coverage. This sets the target for PSF reconstruction for 
a dedicated experiment. In this case the cross-cannel cali- 
bration requirements also becomes stricter by a factor four, 
achievable with current Fourier Spectrometer technology. 
Consider OIII and Nil: the transitions quoted above (88.3 
/im and 205.30 /im, respectively) will dominate the effective 
optical depth to CMB photons because of the pop III yields 
for these elements and because of the transition oscillator 
strength. For these elements, there is an optimal frequency 
range to probe reionization. The upper limit to the redshift of 
reionization sets the reference channel to be at > lOOGHz. On 
the other end, at, say, 350 GHz, the level of contamination 
due to dust in the Galaxy and IR sources becomes important. 
We conclude that the optimal frequency window will be 
between 100 GHz and 350 GHz. This range fits well with the 
frequency coverage of Planck, and, at the same time, is well 
accessible by ground-based experiments. 

In summary, we have presented an algorithm to realistically 
measure the signature of metals in the power spectrum of the 
CMB. We have shown how the different foregrounds can be 
removed and have presented a method to combine the maps 
at different frequencies optimally. The dominant foreground 
is the tSZ and the dominant source of error is the beam error 
and cross-channel calibration. We have forecasted the perfor- 
mance of the algorithm for an experiment with characteristics 
similar to ACT and have shown that one can constrain IGM 
abundances larger than 1% solar at z > 10. Since the method 
probes the metal abundance, and therefore the star forma- 
tion rate, at the redshift corresponding to the frequency of the 
band, it is possible to do tomography of the metallicity enrich- 
ment process. Multifrequency CMB experiments can provide 
a tomographic reconstruction of the re-ionization epoch that 
is nicely complementary to other reionization probes. 
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